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ABSTRACT 

The Chandra monitoring observations of The Antennae (NGC 4038/39) have led to the discovery 
of a variable, luminous, supersoft source (SSS). This source is only detected at energies below 2 keV 
and, in 2002 May, reached count rates comparable to those of the nine ultraluminous X-ray sources 
(ULXs) detected in these galaxies. Spectral fits of the SSS data give acceptable results only for a 
~100-90 eV blackbody spectrum with an intrinsic absorption column of Ah ~ 2 — 3 x lO^^cm"^. 
For a distance of 19 Mpc, the best-fit observed luminosity increases from 1.7x10^^ erg s^^ in 1999 
December to 8.0x10"^^ erg s~^ in 2002 May. The intrinsic, absorption-corrected best-fit luminosity 
reaches 1.4x10*° erg s^^ in 2002 May. The assumption of unbeamed emission would suggest a black hole 
of >100 M0. However, if the emission is blackbody at all times, as suggested by the steep soft spectrum, 
the radiating area would have to vary by a factor of ~10^, inconsistent with gravitational energy release 
from within a few Schwarzschild radii of a black hole. Viable explanations for the observed properties 
of the SSS are provided by anisotropic emission from either an accreting nuclear-burning white dwarf or 
an accreting stellar-mass black hole. 



Subject headings: galaxies: peculiar — galaxies: individual 
— X-ray: binaries sources 



galaxies: interactions — X-rays: galaxies 



1. INTRODUCTION 

At a distance of 19 Mpc {Ho = 75), NGC 4038/39 (The 
Antennae) have long been studied as the nearest exam- 
ple of a galaxy pair undergoing a major merger (Toomre 
& Toomre 1972). In the X-ray band of 0.1-10 keV, the 
first Chandra observation of this system in 1999 Decem- 
ber revealed an extraordinarily rich population of lumi- 
nous point-like sources (Fabbiano et al. 2001). We are now 
in the midst of a year-long Chandra monitoring program 
of The Antennae. The first results of this program, on 
the luminosity and spectral variability of nine ULXs (see 
Fabbiano 1989 and Makishima et al. 2000 for earlier work 
on ULXs), detected with luminosities Lx > lO'^^ergs s~^, 
are reported in Fabbiano et al. (2002). 

Here we report the discovery of a very luminous, vari- 
able, supersoft source in the Antennae galaxies. While 
SSSs with blackbody spectra of ~ 40-100 eV have been 
detected in several galaxies (e.g., in Local Group galaxies: 
Greiner 1996; M81: Swartz et al. 2002; MlOl: Pence et 
al. 2001), their typical luminosities do not exceed 10^^ erg 
and are mostly lower than 10"^^ erg s~^. These sources 
are believed to be the result of nuclear burning on the sur- 
face of accreting white dwarfs in binary systems (van den 
Heuvel et al. 1992; see also Kahabka & van den Heuvel 
1997). However, the source of interest here, CXOANT 
J120151.6-185231.9 (= Source 13 of the Zezas et al. 2002a 
list), may reach luminosities in excess of 10*° erg s~^. 

Source 13 was first detected in 1999 December with a 
count rate '^4 times below that of the ULX range. It 

^ http:/ /asc. harvard.edu/cal/Acis/Cal_prods/qeDeg/ 

^ http:/ /www. astro. psu.edu/users/chartas/xcontdir/xcont. html 



is indicated by a circle in fig. 1. Our monitoring of The 
Antennae shows that this source reached a count rate com- 
parable to those of the ULXs in 2002 May, while keeping 
an unusual, very soft emission. In this paper, we present 
both the light curve of this source and a spectral study of 
its emission. 

2. OBSERVATIONS AND ANALYSIS 

Table 1 summarizes the log of the four Chandra ACIS-S3 
(Weisskopf et al. 2000) observations of The Antennae dis- 
cussed in the present paper, and lists the observing times 
after screening for background flares. Details of the data 
analysis are given in Fabbiano et al. (2003). This analy- 
sis includes astrometric correction of the 2001 December, 
2002 April, and 2002 May observations to the 1999 De- 
cember coordinates, and source detection in four spectral 
bands following the prescriptions by Zezas et al. (2002a): 
Full band (0.3-7.0) keV; soft (0.3-1.0) keV; medium (1.0- 
2.5) keV; and hard (2.5-7.0) keV. The data were corrected 
for spatial and spectral variations of the ACIS-S3 response, 
including the time- variable ACIS-S3 effective area^ ^. This 
correction results in a factor of 1.8 increase of the count 
rate in the last three observations, when compared with 
the 1999 December data. The back-illuminated ACIS-S3 
CCD is not affected by the energy response degradation 
experienced by the front-illuminated CCDs; therefore the 
CTI correction is not relevant here. 

Figure 2 shows the light-curve of Source 13 in the 0.3- 
7.0 keV band. We observe a ~8-fold increase of the cor- 
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rected count rate from 1999 December to 2002 May. While 
detected at sub-ULX count rates in the 1999 December 
data (Zezas et aL 2002a), in May 2002 the count rate is in 
the range of the ULX count rate in The Antennae (corre- 
sponding to Lx > 10"^" erg s^^ , for a 5 keV bremsstrahhmg 
spectrum), but the spectrum is much softer than the spec- 
tra of the other ULXs. This source remains undetected in 
the hard band, and only the color {S — M/S + M) can be 
derived. This color is 0.88±0.04, significantly softer than 
the typical ULX colors that all lie in the 0.2 to —0.4 range 
(Fabbiano et al. 2003). 

We used XSPEC for the spectral analysis (see Zezas 
et al. 2002a). Source counts were extracted from cir- 
cles of 3 pixel (1.5") radius, with backgromid estimated 
from surrounding annuli. The data from each observa- 
tion were fitted with a variety of models, available in 
the XSPEC library, including blackbody, disk-blackbody 
(used in the ASCA studies of ULXs, sec Makishima et 
al. 2000), brcmsstrahlung, Raymond-Smith thin- plasma 
model (RS), and power-law. A variable absorption col- 
umn was included in each model. Spectral analysis was 
performed on data at energies 0.3-7.0 keV for all the data 
sets, binned so that there were at least 25 counts in each 
bin. Typically, extremely few spectral counts were de- 
tected above 2 keV in the 2002 May data, and above 1 keV 
in the earlier observations, when the source was fainter. 
This binning results in four bins for the 1999 December 
observation, when the source was faintest. To verify that 
the determination of A^h was not affected by calibration 
uncertainties, we also (a) ran our spectral analysis exclud- 
ing data below 0.6 keV (to avoid the Oxygen edge) and 
(b) added an edge to the models. In all cases, we obtained 
consistent results. 

The RS model is a bad fit for all the data sets. Black- 
body, disk-blackbody and bremsstrahlung all give similar 
quality fits. Best-fit kT is in the range of ^^150-90 cV for 
the blackbody model, 110-220 eV for the disk-blackbody, 
and 120-370 eV for the bremsstrahlung. The power-law 
model is acceptable for the first and third datasets, but 
gives a bad fit for the higher statistics 2002 May data; 
it is also marginal for the 2001 December data. In gen- 
eral power-law indeces tend to be very large, refiecting 
the very soft spectra. These results, including the best fit 
and the degrees of freedom (number of bins minus fit 
parameters) are summarized in Table 2, where the errors 
are Irr for 1 interesting parameter. To verify that diff'er- 
ent binning did not affect our results in our lowest count 
rate dataset (1999 December), we also analyzed these data 
with 15 counts per bin, as in Zezas et al. (2002a), obtaining 
consistent results. 

The best-constrained results are from the 2002 May 
observation, when the source was most luminous: the 
temperature then was kT ~ 90 eV for the blackbody 
model, 110 cV and 120 eV for the disk-blackbody and the 
bremsstrahlung models, respectively. The fits also suggest 
a large absorption column Ah ^ 3 — 5 x lO^^cm"^, de- 
pending on the model. For the rest of this paper we will 
assume that the emission is optically thick, and for simplic- 
ity use the blackbody model. Use of the disk-blackbody 
model would not change our conclusions. Figure 3 shows 
the data, best-fit blackbody models, and residuals. Fig- 
ure 4 shows the la 2-parameter Nu-kT confidence con- 



tours. For the 2002 May data we also show the 99% con- 
fidence contour. Fitting a two-component model (black- 
body or disk-blackbody plus power-law) did not produce 
a significant improvent of the fit. 

The observed best-fit source luminosity (0.1 2 keV) 
varies from 1.7x10'^* erg in 1999 December to 

8x10^* erg s'^ in 2002 May (Table 3). These values can 
be considered to be lower limits to the intrinsic source lu- 
minosity, since no extinction correction has been applied. 
The intrinsic (i.e., emitted) best-fit luminosity is signifi- 
cantly larger in 2002 May, because of the large Ah required 
by the fit, and reaches 1.4x10^° erg s~^. In 2002 May, we 
estimate a minimum intrinsic luminosity of 4 x 10^^ erg 
by calculating the flux for the la lower limit on both 
kT and Ah- If the emission of the SSS is due to nuclear 
burning on a white dwarf (WD) surface (van den Heuvel 
et al. 1992), a hot WD atmosphere may be more physical 
than a simple blackbody. However, based on the results 
of Swartz et al. (2002), we estimate that the adoption of a 
model WD atmosphere will result in less than a factor of 
2 change in the estimated luminosity. 

3. DISCUSSION 

What is this supersoft, luminous, variable source? The 
large absorption column and the lack of an obvious iden- 
tification in the HST WFPC2 data (Zezas et al. 2002b) 
exclude a foreground object. Given the low density of 
supersoft AGNs (Puchnarewicz 1998), the likelihood that 
this source could be a background object is low. If this 
source belongs to The Antennae, it is the most luminous 
galaxian SSS ever detected. Our observations suggest that 
the source has a thermal spectrum. Although the uncer- 
tainties are large, the best-fit temperature does not seem to 
increase with luminosity. Taken at face value, this behav- 
ior is incompatible with isotropic emission from a constant 
radiating area. In the following analysis, we investigate 
these trends systematically. 

King fc Puchnarewicz (2002) show that blackbody emis- 
sion from a region of size r times the Schwarzschild radius 
of a mass M obeys the relations 
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where Lsph is the inferred isotropic luminosity of the 
blackbody, L the true source luminosity, b the beaming 
factor that accounts for eventual non-spherical emission, 
/ = L/L^^dd the Eddington factor, p ~ 1 a measure of the 
geometrical deviation from a spherical photosphere, and 
Tioo the source temperature T in units of 100 eV. 

Dividing (1) by (2) gives Im/b in terms of igph, where 
m = M/Mq. Dividing (1) by the square of cq. (2) gives 
pr'^m? /b in terms of LgphT^*, and we can deduce l/mpr^ 
in terms of T^. We calculate these quantities (Hsted 
in Table 3) from the results of Table 2 for the black- 
body spectrum. It is easy to show from Table 1 that 
[Im/b) oc {jtrr'^m'^ /bY/^, so that 

r ^ ;3/4p-l/2^-l/4^ 

Thus at most one of the three quantities b,l,pr'^ can be 
constant. We consider three cases. 
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(i) Near-isotropic emission from an interm,ediate-mass 
black hole (6 ~ 1 ~ p). Given the peak Imninosity of 
the SSS, this case requires a black- hole mass >100 M0(cf 
Miller et al. 2003) However, the large increase (a factor 
^1000) in radiating area is very hard to reconcile with a 
simple picture in which a black hole accretes from an accre- 
tion disk and the blackbody emission comes from a region 
of a few Schwarzschild radii. Colimm 3 of Table 3 shows 
that I must increase by a factor 81 between the 1st and 
4th observations; for this type of geometry we clearly have 
p ~ 1, so the radius factor r must increase by a factor ~25. 
Since the blackbody emission comprises most of the puta- 
tive accretion luminosity, it must come from deep within 
the potential well, i.e., r < a few. So, if we (generously) 
set r ~ 3 in Observation 4, the emission must be confined 
to an implausibly small region of only 0.1 Schwarzschild 
radii in Observation 1. 

Although the backbody fits point to the exclusion of 
an massive black hole, this conclusion is not iron-clad. It 
hinges on the assumption that the 1999 December emis- 
sion is indeed a blackbody. As discussed earlier, this is the 
least well constrained dataset, so that we cannot exclude 
that the spectrum may follow a power-law distribution. 
Although the power-law F ~ 4 we obtain for the 1999 De- 
cember data is extremely steep, in excess of typical black- 
hole binary power-laws, the uncertainties are large (see 
also Zezas et al 2002a). Power-law components have been 
seen to dominate the emission in ULXs in low state (see, 
e.g. La Parola et al. 2001; Kubota et al. 2001; Fabbiano 
et al 2003). If the first observation is discarded from the 
blackbody analysis of Table 2, the required radius increase 
between observations 2 and 4 is reduced to a factor ~ 2. If 
we only consider observations 3 and 4, we obtain a radius 
increase of a factor of ~ 4. It remains to be demonstrated 
how this smaller increase would fit into an intermediate- 
mass black hole picture. 

(a) Varying beaming at constant luminosity. An oppo- 
site extreme from the near-isotropic case is I ^ constant, in 
which case b decreases by a factor ~139 between Observa- 
tions 1 and 4. The last column of Table 1 then shows that p 
must simultaneously increase by a factor ^10. Physically, 
holding I constant while other quantities vary widely is 
plausible only in one situation, namely that the source is 
radiating constantly at the Eddington limit (^ = 1) while 
the accretion rate may change. In this case the range of b 
is m — m/185. The l/m,pr'^ column of Table 1 then shows 
that pr^ = 2.7 x lO^rn"^ in Observation 4, giving the ra- 
diating object a radius of i? = 4.8 x lO^m^/^p"^/^ cm. For 
m ~ 1, p ~ 1 this radius is suggestively close to the radius 
of a white dwarf. This is not surprising, as the inferred 
temperature and luminosity are now typical of supersoft 
X-ray binaries, which are thought to be powered by nu- 
clear burning of matter accreting on to a white dwarf. 

If the SSS is an accreting WD. our results would be 
consistent with an increasingly super-Eddington accretion 
flow (M > MEdd ~ 10"'' M0yr~^ for steady nuclear burn- 
ing on a WD surface). As M increases beyond MEdd, the 
flow geometry apparently changes gradually from a thin 
disk (1999 December), with burning on a narrow, isotrop- 
ically emitting {b ~ 1) equatorial band {p ^ 0.1) on the 
white dwarf, to a thick disk with burning on most of the 
WD surface (p 1), but now with a strongly anisotropic 



radiation pattern (6 ~ 10^^) in our direction. The most 
likely cause of this anisotropy is warping of the accretion 
disk. We note that disks are known to warp in supersoft 
X-ray binaries (Southwell et al. 1997), essentially because 
of their proximity to the Eddington limit. At still higher 
M nuclear burning on white dwarfs drives either envelope 
expansion where the source swells up to red-giant dimen- 
sions, or a vigorous wind outflow, or both (cf Hachisu et 
al., 1996). In all cases the system is likely to be extin- 
guished as an X-ray source. 

(Hi) Mildly anisotropic emission (b > 0.1^ from, a 
stellar-mass black hole. At first sight the pliotospheric 
radii of ^10^ -10^ cm we deduce from our observations do 
not appear any more natural for a stellar-mass black hole 
than for an intermediate-mass one as in (i) above. How- 
ever, Mukai et al. (2002) have pointed out that accretion 
at rates comparable to Eddington must lead to outflow, 
and have shown that the electron scattering opacity of 
the resulting wind does imply supersoft emission with a 
photospheric size of this order. The MlOl source studied 
by Mukai et al. (2002) has a supersoft luminosity of order 
2q39 ^-|.g g-i hence, does not require anisotropic emis- 
sion for a black hole mass >10 Mq. However, their analysis 
is easily extended to the case that an Eddington-limited 
source blows out a wind conflned to a double cone of to- 
tal solid angle 47r5 about the black- hole axis. Since this 
wind represents the path of lowest optical depth through 
the accretion flow, the radiation will escape this way also, 
implying p ~ 6. We follow Mukai et al. (2002) in as- 
suming a constant velocity for the outflowing material, 
since this material is likely to achieve escape velocity and 
coast thereafter. We neglect any emission from this wind 
and compute the Thomson optical depth t by integrating 
the electron density from radius R to infinity. Since 
/ Nedr ~ Mout /47r6wi?mH we find a photospheric radius 
3x10^ . 

i?ph = — Mi9 cm, (4) 

bug 

where Wg is v in units of 10^ cm s~^ and Mig is the outflow 
rate in units of 10^^ g s^^, the Eddington accretion rate 
for a 10 Mq black hole. Clearly, we must choose a mass of 
this order for consistency. The luminosity ~10^^ erg s~^ 
of Observation 1 is then definitely sub-Eddington, so pre- 
sumably there is very little outflow and we see down to the 
inner accretion disk directly in this observation. Observa- 
tions 2, 3, and 4 are all close to ^Edd. i-c, 1 = 1. Since 
b = p and the true blackbody radius i? = 3 x lO^rm cm, 
we can read off the size of the photosphere directly from 
the last column of the Table, and l/b also. This gives an 
inner disk radius of order 2 x 10* cm for Observation 1, 
with i?ph = 2 X 10^ cm, 5 x 10^ cm, and 5 x 10^ cm for 
Observations 2, 3, and 4, respectively. Simultaneously l/b 
increases from ^^0.1 to 3.1, 0.6, and 11. This is consistent 
with I = 0.1, 1, 1, and 1 over the four observations and b 
eventually decreasing to a value of ~0.10, presumably as 
Mace rises above the Eddington rate. Self-consistently, the 
assumed vg is above the escape value for the values of Rph- 
Of course, b would be larger still if we took a BH mass of 
15 Mq rather than 10 Mq, as observed in GRO J1915-I-105 
(Greiner et al. 2001). 

The above values are consistent with the suggestion by 
King et al. (2001) and King (2002) that ULXs are actu- 
ally X-ray binaries involving stellar-mass black holes, but 
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with mildly anisotropic radiation patterns (b ~ 0.1) re- 
sulting from accretion at close to the Eddington rate onto 
the black hole. The above scenario can arise in two cases 

(i) thermal-time-scale mass transfer when the companion 
star in a high mass X ray binary fills its Roche lobe, and 

(ii) bright outbursts of soft X-ray transients (King, 2002). 
The first case is dominant in most galaxies, although the 
second must account for ULXs in ellipticals. Confirma- 
tion that this is a reasonable explanation for ULXs as a 
class comes from Grimm, Gilfanov, & Sunyacv (2002), who 
show that — when normalized by the star formation rate — 
ULXs form a natural extension to the luminosity function 
of high-mass X-ray binaries in nearby galaxies. 

4. CONCLUSIONS 

We have discovered a variable SSS {kT ~ 90 keV) in The 
Antennae, Source CXOANT J120151.6-185231.9, which 
reached a peak intrinsic luminosity of 1.4 x lO"**^ erg s~^ in 
2002 May. 

Near-isotropic emission from an intermediate-mass 
black hole accreting from a disk would be incompatible 
with our observations of this source in the most likely case 
of soft thermal emission, as the radiating area would have 
to increase by more than a factor 1000 over the four ob- 
servations. There remains a less likely possibility that the 
1999 December emission may be due to a low-intensity 
power-law dominated state (e.g. Kubota et al. 2001), in 
which case the required area increase is reduced to a less 
demanding factor 10. 

A possible solution is a white dwarf with M ~ 1 Mq, 



accreting at the Eddington limit (1 = 1) and with a vari- 
able beaming factor. This explanation has the advantage 
of giving a natural scale for the deduced photospheric ra- 
dius, but does require extreme beaming (up to 6 ~ 10~^) 

A second possible solution involves outfiow from a 
stellar-mass black hole, accreting near the Eddington limit 
(Mukai et al. 2002). A consistent explanation of our ob- 
servations results if this hypothesis is combined with the 
suggestion by King et al. (2001) that ULXs are actually 
X-ray binaries involving stellar-mass black holes, but with 
mildly anisotropic radiation patterns (6 ~ 0.1). 

The stellar-mass black hole solution is the more conser- 
vative choice in the present case, as it docs not require 
as extreme an anisotropy as a white dwarf. However the 
latter may remain a realistic candidate for slightly less lu- 
minous supersoft ULXs. 
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Table 1 
Observations 



OBSID 


Date 


Net Exp. 






(Ks) 


315 


1999-12-01 


75 


3040 


2001-12-29 


64 


3043 


2002-04-18 


61 


3042 


2002-05-31 


67 



Table 2 



Spectral Fits 



OBSID 




B. Body 






Disk-BB 






Brcmss. 






RS 






PO 






kT 


Nh 




kTi„ 


Nh 




kT 


Nh 




kT 


Nh 


x' 


r 


Nh 


x' 




keV 




d.o.f. 


keV 


10^^ cm~^ 


d.o.f. 


keV 


10^^ cm~^ 


d.o.f. 


keV 


10^^ cm~^ 


d.o.f. 




10== cm~ = 


d.o.f. 


315 


(, 1C + 0.03 
"•■'^O_0.04 


< 0.12 


1.3 




< 0.15 


0.43 


f, ,7+0.27 


0.04(< 0.2) 


0.3 


0.22 


< 0.04 


13.4 


4 74+5.3 
'^-1.3 


18+°" 


0.23 






1 




1 




1 






1 




1 


3040 


"•^"-0.024 




4.0 
4 


11+004 


0.25iS:f, 


3.8 
4 


Q 14+0.05 


n 00+0.32 


3.6 
4 


u.uy_Q Q3 


0.63tg:l? 


11.4 
4 


9.6(> 7.8) 




6.18 
4 


3043 


f, -,0+0.02 


0.025(< 0.22) 


2.2 


14+°°^ 


0.09(< 0.29) 


1.91 






1.64 


0.20 


<0.1 


18.5 


2.641=,;° 


46+° °* 


2.17 






3 




3 


3 






3 


3 


3042 


09+° °^ 




21.7 
10 


11+°°^ 




22.4 
10 


12+°°^ 

"•^•^-0.03 


0.49lS:l? 


23.5 
10 
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Table 3 

Blackbody Luminosities and Derived Parameters 



OBSID 
(1) 


Lx. 0.1-2 kcV 

10^" org s"^ 
(2) 


T omitted 
^X, 0.1-2 koV 
lO-'** org S"^ 
(3) 


Im/b 
(4) 


/ / mpr^ 
(5) 




(6) 


315 


0.17 


0.17 


1.33 


2.8 X 10" 


-6 


4.8 X 10° 


3040 


0.43 


3.9 


30.5 


5.6 X 10" 


-7 


5.4 X lO'' 


3043 


0.53 


0.8 


6.2 


1.6 X 10" 


-6 


3.9 X lO" 


3042 


0.80 


13.8 


108 


3.6 X 10" 


-7 


2.9 X 10* 




■ 





Fig. 1. — The Dec 1999 image of The Antennae (Fabbiano et al. 2001), with the super soft ULXs discussed in this paper identified by a 
circle and source number from Zezas et al. (2002a). Soft emission is red, hard emission is blue. Luminous sources appear white because they 
emit in the entire spectral band. 
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1 2 - 20 ks 



Dec 99 



Dec 01 Apr 02 



Src 13 



May 02 



Fig. 2. — The light curve of CXOANT J120151. 6-185231. 9. The starting date of each observation is given in Table 1; the data are binned 
in 10 ks time bins. 




Fig. 3. — Observed spectral data (see text) and best-fit blackbody spectra, with the fit residuals. Black: 1999 December; green: 2001 
December; blue: 2002 April; red: 2002 May. We show only the 0.3 - 2 keV range, where the greatest majority of counts were detected. In all 
cases the last bin used for the fitting extends to 7 keV. 
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Fig. 4. — Nn-kT confidence contours. Sofid contours are at la for two interesting parameters. For the 2002 May observations, the dashed 
contour is at 99% for two interesting parameters. 



